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Abstract 

The photon structure function is a useful testing ground for QCD. It is perturba- 
tively computable apart from a contribution from what is usually called the hadronic 
Qh! component of the photon. There have been many proposals for this nonperturbative 



part of the real photon structure function. By studying moments of the virtual pho- 
ton structure function, we explore the extent to which these proposed nonperturbative 
contributions can be identified experimentally. 

1 Introduction 

With the discovery of asymptotic freedom, QCD has become the accepted theory of hadrons. 
However, at the energy scale of the formation of hadrons (about a few GeV or less), per- 
turbative QCD (PQCD) can no longer be applied alone, and nonperturbative QCD (NP) 
becomes important and may even rule. In this regime, we have very limited tools we can 
use to understand what is going on. 

There is a fundamental dimensional constant giving a qualitative citerion of separating 
these two energy regions, = A. We use the modified minimal subtraction renormaliza- 
tion scheme, MS, throughout this paper. We can hope for some progress in understanding 



the interplay of nonperturbative and perturbative QCD in those processes where both play 
a role. It may even be important to understand the relation of perturbative and nonpertur- 
bative contributions in order to obtain A from experiment. 

The photon structure function in e^e~ — > e^e~ + hadrons has been a useful tool to 
study QCD. It has been calculated perturbatively to nonleading order, which is necessary 
to obtain a meaningful comparison with QCD predictions for other processes While 
valid asymptotically, the PQCD prediction is modified at attainable by the presence of 
a NP part, usually called the "hadronic" component of the photon. (This is an unfortunate 
usage and we will refer in the following to the PQCD and NP contributions to the photon 
structure function.) Our point of view in this paper is that the photon structure function 
is valuable precisely because it contains both PQCD and NP contributions at not too large 
Q^. This may make it a unique process in which we can study both of these as functions of 
the kinematic variables. 

The most information can be extracted from this process by studying the photon structure 
function for a virtual photon target whose invariant mass, "— p2"^ f^^^^ vary over a range from 
zero to some appreciable <^Q^. 

Here we exploit the moments of the photon structure function from the calculation of 
Uematsu and Walsh for the perturbative part of the struncture function 0. (In particular, 
we use the normalization and conventions of that paper.) We use the suggestions of various 
authors for the nonpurturbative part [Q, [Q. This "hadronic" part to the function has 
been obtained by these authors by fitting experimental data from other processes. Typically, 
one finds the pion structure function from a semihadronic process such as vr + p — > /i+yU^ + 
anything. The static quark model and vector dominace are then used to obtain a conjectured 
NP component to the photon structure function. 

The basic idea of this paper is to examine the behavior of the total photon structure 
function (including both PQCD and NP parts) as a function of and P^, with the aim 
mentioned above of finding out how to distinguish the two experimentally. First we describe 
models for the NP contribution-including the dependence often left out of consideration. 
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We then discuss the real and virtual photon structure functions. We then present our results 
for the behavior of the two components of the structure function. 

We are interested in the kinematic region 2 — 4 GeV^ < < 50 GeV^ and < < 1 
GeV^. We would like to stimulate experimental studies of this range of variables. 

Throughout this paper, we set the quark mass to zero and the number of quark flavors 
to four. While interesting and relevant to a comparison with experiment at high Q^, a 
consideration of quark mass effects here would distract from our main point. 

2 Models for the Nonperturbative Component 

We think that it is very much an open question how the nonperturbative component of the 
photon structure function is to be treated. For the present, we describe the conventional 
view. We return to the question later in the paper. According to this conventional view, the 
overall real photon structure function is a sum of two parts, 

F^{x, Q') = FIpqcd{^, Q') + F^^pix, Q'). (1) 

It is customary to write the nonperturbative contribution to the photon structure function 
^2NP ^ coherent sum of light vector meson states. From this vector meson dominance 
model we can then identify the NP part of the photon structure function 0], approximated 
by the p meson contribution, 

FlM^,Q')c:.{^)Ff{x,Q'). (2) 

Where 7^ = /^/tt ^ 2.2 and F2 can be obtained from the pion structure function using the 
quark model 

Ff{x,Q') = Ff{x,Q') = + ^f,/^-) = ^xf^/^-. (3) 

Here we denote xfuf-n- as xu{x, Q^) for convenience and similarly for the others. We also 
put u{x, Q"^) = d{x, Q^) in the pion. This assumes that we are at low enough so that we 
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can neglect sea quarks. Then equation becomes 

FImpM') = (^)Ixu{x,Q% (4) 

From now on we confine ourselves to < 50 GeV^ where we expect that our neglect 
of sea quarks will be adequate. From definitions of the non-singlet and the singlet quark 
distribution functions we can write a quark structure function in terms of the non-singlet 
and singlet parts of the structure function [^]. 

xq,{x, Q2) = xqf'ix, Q^) + Q^), (5) 

where = J^iiQi + Qi) ^^nd i runs over all available flavors, here we choose the number of 
flavors to four, / = 4. The final form of F^^^fp of the photon structure function becomes 

^Inp = (^)^Nr(^,Q') + ^xq'{x,Q% (6) 

We have xq['^ = xu'^'^ for the pion in this example. The moments of the nonperturbative 
part of the structure function are then 

F^^pin, Q') = f dxx^-'Fl^p{x, Q'). (7) 

JO 

Our procedure will be to obtain xq'^^ and xq'^ at all through numerical integration of the 
Altarelli-Parisi equations. So we can calculate M^j^p{Q'^) at any Q^, given its starting value. 
We now need to determine the starting value of the p meson or pion structure function at 
some low Q^. 

We use three different fits which were obtained by three groups Newman et. al. [Q, NA3 
collaboration at CERN and Castorina et al. |p. All these fits have been obtained at 
(5^ — 2 — 4 GeV^. We will show all the solutions to the Altarelli-Parisi equations for the 
corresponding NP contribtutions to the structure function. We emphasize that we have no 
convictions as to which of the resulting NP parts is the correct one. Our aim later in this 
paper will be to see to what extent they can be distinguished experimentally for > 0. 

We will include three examples of Non-Singlet (NS), Singlet (S), Gluon (G) parts of the 
structure function at Qq. 
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The general form of the parametrizations of the hadronic functions can be written as 
xu — Cx^-ll — xY since the three fits are given as 

' (A) Newman et al. : xu{x, Ql) = 0.52(1 — x) 

< (B) NA3 coUabaration : xu{x, Ql) = 0.57a;°-^(l - x) (8) 
, (C) Castorina et al. : xu{x, Ql) = 0.766x^-^{l - xf-^ 
which yield the nonperturbative structure functions considering only the dominant p meson 

0.13 a (1-x) (A) 

0.14 a x^-\l - x) (B) (9) 
0.19 a x^-^{l - x f^ (C). 
The first of these is slightly at odds with our neglect of sea quarks, since it goes to a 
finite limit as x — > 0. Since our aim is only to compare different assumptions about the p 
meson structure function for high moments, we will use it anyway. 

Let us consider the ratios of the momenta carried by quarks and gluons. Since u — d 
inside the vr" at low Q^, we can only consider the u for quarks. We define the ratios for the 
quarks and gluons as 

rl 

ry. = lo dxxf^i^- = / dxxu{x, Q^) 

Jo 

rl 

Tg^ J^dxxfgi^-^ I dxxg{x,Q'^). (10) 

t/ 

We can calculate using the three fits 



r^^C C dxx^il - xf = S(a + 1, 6 + 1), 
Jo 



(11) 



where B is the Euler Beta function with Re (a+ 1) > 0, Re ip+l) > 0. Therefore the values 
of for each fit are 



'0.52fi(l,2) =0.26 (A) 

ru = < 0.57 5(1.4,2) = 0.17 (B) 

.0.766 5(1.5,1.6) =0.28 (C). 
Prom the QCD momentum sum rule we have rg — 1 — 2ru so that 

r 0.48 (A) 

= 1 -2 C5(a + 1,6+ 1) = i 0.66 {B) 

io.44 (C). 



(12) 



(13) 



5 



We can determine directly from ru with the assumption of no sea quark contribution 
to the total momentum at the low initial value of regardless of the form of the gluon 
structure function. The parametrization of xg{x, Qq) is up to us to determine. The simplest 
form of the gluon structure function at is 

xgix,Ql) = C^il-xf\ (14) 

We can find the value of Ci from the QCD momentum sum rule after making an assumption 
of the value of C2 {C2 > 0) which could in principle be determined from experiment. We 
choose C2 = 2 here. Using the momentum sum rule again we obtain 

= JTTC^) =^-^C B{a + l,b+l). (15) 

This equation yields 

r 0.624 (A) 

C2 = 2, Ci = \ 1.584 (B) (16) 
[ IMA (C). 

Let us look at the values of Vu for the three examples. In example (A) 26% of the total 
momentum is carried by the u quark and thus 48% of it is carried by gluons in p (or 7r)-meson. 
Similar agument applies to other examples although, with the model that we suggest, the 
example (B) seems to have too big gluon momentum (66%) compared to quarks (34%). The 
authors of reference (C) chose xg{x, Qq) = 1.408(1 — x)^. Then = 0.35 which seems small 
compared to our result 0.44. However, one needs to keep in mind that the model for the 
gluon structure function here is only the simplest one and all three fits to the experimental 
data are perfectly acceptable. 

This analysis is independent of the number of vector mesons. However, in the nonpertur- 
bative part of the structure function we need to sum over all the vector mesons coherently. 
For example, we will have about 60% enhancement in the nonperturbative part when we 
include the incoherent sum of all vector mesons u and in addition to the dominant p [0. 
In this paper we confine ourselves only to the p-term for simplicity. 

In the numerical calculation, we evolve the equations from Qq = 3 GeV^ to = 45 
GeV^ with A = 0.2 GeV. The evolution of the nonsinglet function can be obtained in a 
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straight forward manner using the above hadronic functions. In coupled equations of singlet 
and gluonic functions we assume the initial gluon distribtuion in the form of equation ( |T4|) 
with no initial sea quark distribution. We use the normalization factor a ln((5^/A ) and we 
cut at very low x [x < 0.025) and high x [x > 0.975) to avoid numerical problems at the 
end points. 

We show the moments of the NP part of the structure function, F^^^pin^ Q^), in Figure |I| 
with the normalization factor, ot ln((5^/A ). The qualitative behavior of NP of three fits are 
same. They decrease fast at low while the rate of the decrement becomes small at high 
. This means that we do not expect very small NP contribution at a given moment to the 
structure function even at high . Therefore we expect the NP contribution to the total 
structure function persists even at high as we will see in the next section. We will come 
back to this point in the real and the virtual cases. 



3 The Real Photon Structure Function 

The PQCD moments of the real photon structure function can be written as 

F^,PQCDin,Q')= dxx--'F^{x,Q') = ( a„ + 6„) 

Q2 ^ q2 

~ — (a„ ln=2 + ctn In ln=2 + &n)- (17) 
A A 

Where the expressions and the values of a„, 6„, and a„ are given in reference [Q]. The 
factor 1/e^ in the right hand side of the equation is from the difference in normalizations of 
structure functions between reference [0 on which our calculation is based and reference [|1|. 
The approximation comes from omitting the higher orders in 
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Po In^ O 
A 



lnln(gVA' 

?2 , /^o ,-r2s 



^Po ln(Q2/A' 

A is given in the MS scheme and g'^{Q'^) = ^|jg(<5^) is the coupling constant with two- loop 



order term 



tm 1 Pi In ln(gVA' 



167r2 ln(Q2/A') Pi In^iQ^A'" 



(19) 



Let us go back to the equation (0); then we have 

an = ^CLn- (20) 

Po 

We spht the total structure function into two parts as before. We use the resuhs in 
for the perturbative QCD part of the structure funciton and the previous three examples 
for the NP part. We can recover Bardeen and Buras's moments from Uematsu and Walsh's 
ones formally by putting = A . 

We write the moments of the real photon structure function in equation (||) as 

F^{n,Q') = F^^pQCDin,Q')+FlMn,Q'). (21) 

We obtain the total real photon structure function in each moment by summing the 
nonperturbative and the perturbative contribution for the corresponding expressions. As we 
saw in Figure 0, the NP parts of the structure functions decrease relatively fast at low 
while the rate slows down as increases. On the other hand PQCD increases fast at low 

and slowly at high as we see in Figure 0. This is a reasonable result since we expect 
that NP contribution becomes larger than that of PQCD at a low Q^. If we add the two 
contributions, the total structure functions, F2tot'''{n,Q'^), for n = 4 increase very slowly at 
high after normalizing by the factor a \n{Q'^/A^). Here the dashed curve indicates the 
PQCD part and A, B, and C indicate the total structure functions for the three examples. 
All the moments (for n > 4) increase rather slowly after > 15 GeV^. We choose different 
scales for Figure ^(a) and Figure |^(b) to distinguish explicitly the curves for those examples. 
We also notice the faster decrease of the NP part compared to PQCD at higher moments. 
The total real structure functions for all three fits are more or less the same. We found that 
the NP contribution to the total structure function is about 10% around = 3 GeV^ when 
n = 4 and slightly less when the moments become higher and gets larger. 

4 Virtual Photon Structure Function 

Let us now move on to a virtual target photon. Consider a deep inelastic scattering involving 
two photons with momenta q for the probe photon and p for the target photon. We have 
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'Q^ and = —P^ as the probe photon mass and the target photon mass respectively. 



We take the kinematic region of the momenta of the photons as follows. 



-r2 



A < < 



(22) 



The moments of the structure function can be calculated by solving the renormalization 
group equation for the Wilson coefficient in the operator product expansion of photon-photon 
scattering to leading order in a and to next to leading order in as- 



^2,PQCD 



Jo 



167r2 2/3o 



y- pn_ 



1 



+ Af/2/3o^2(g2 

tin) 



A7/2A,+1 



(23) 



where P", A", P", C" and A" are given in references 



and the index i runs over 



+, — , NS. Notations, +, — , NS are from the eigenvalues, A!^, A^_5, of the one-loop hadronic 
anomalous dimension matrix. tiQ"^) is the coupling constant to two loops. Now, equation 



(p3|) is valid forn = 2 since the singular behavior of the coefficient A" , which causes a problem 
in the real photon structure function, is cancelled by the factor, 1 — ((7^((5^)/(7^(P^))'^"^^'^" 
when A" goes to zero. Since the expression (p3D is independent of any renormalization 
scheme, to change the equation from one scheme to another all we have to do is to write g, 
A^, Pf, and A in the new scheme. 
We now write 



F^ix, Q\ P') = FlpQani^, Q\ P') + P^^x, g^ P'^ 



with the corresponding moment functions 



P2^(n, Q\ P') = F^pQc^in, Q\ P') + P^^n, Q\ P'). 



(24) 



(25) 



We continue to follow the convention and the simple vector meson dominance model for the 



nonperturbative part of the structure function by writing its dependence on as 

F^M-,Q\P') = ^^^^^. (26) 

Where -F^tvp ^^^^ photon nonperturbative part which we have been discussing up to 

this point. In the previous section we simply had = 0. This is the only consistent way 
of including the P^ dependence once we adopt the notion that the NP part of the photon 
structure function can be approximated by the simple vector dominance model. 

As n increases, the moments are sharply suppressed in exactly the same way as for the 
real photon case. Note that the n = 2 case is possible for both the nonperturbative and 
perturbative parts in the virtual photon structure function while only the nonperturbative 
part appears in the real structure function case. 

Since the PQCD structure function of the real photon can be obtained from the virtual 
photon structure function formally by putting = A , it makes sense to replace P^ in our 
expressions by + A . Then we return to the real photon structure function by setting 
P^ = rather than A . This modification does not affect our calculation much since we need 

2 

P^ = and P^ ^ A . Henceforth we write the PQCD piece with this change. 

In Figure ^ we show the total virtual photon structure function, P^jo^l^^ Q^y -P^); versus 
the probe photon mass from PQCD and NP after normalizing the structure function 

2 

by a ln((5^/(P^ + A )). These display the strong nonperturbative dependence of the total 
structure function at low {Q"^ < 10 GeV^) and at low moments. Perhaps the most 
interesting behavior of the function occurs in the regions 0.4 GeV^ < P^ < 0.6 GeV^ and 
Q > 20 GeV^. As the moments increase we observe similar behavior to that in the real 
photon case. 

It is useful to look at the n = 2 moment in the virtual structure function since there 
is no such moment in the real case. In fact we found interesting behavior of the structure 
functions in our three examples. When P^ = 0.2 GeV^, we have qualitatively same behavior 
for all three fits while the fit (A) behaves very differently when P^ = 0.4 GeV^. However, 
one should keep in mind that we have used only p meson in the analysis. 
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From these plots we can clearly see how the NP contribution behaves, given different 
assumptions about its form at = 0. 

We close with a comment on the PQCD transition to the real photon case. The virtual 
photon structure function can be reduced to the real structure function by taking the target 
photon mass equal to zero (i.e., = A in equation (P3|)). Then g^{P^ = A ) goes to the 
infinity where quark confinement start occurring. The equation p3|) becomes 



Jo 



1677^ 



(27) 



167r2 2(3o 

Now this expression is valid only for n > 2 since we cannot calculate the moment for n = 2 
case in the real structure function due to the existence of unknown hadronic matrix elements. 
If we compare two equations (p!7| ) and (^) in the MS scheme (A and g^j^ are defined in the 
equation (p!9|)), 



«M67r2 



167r2 2po 



167r^ 



X2/2(3og\Q': 



+ E^r + E5r + q 



(2^ 



From the above equation we can easily identify 



bn 



1 

Wo 

1 

Wo 



+ \-/2(3otm 



(29) 



Therefore we obtain all coefficients in Bardeen and Buras' form from those in Uematsu and 
Walsh. 



5 Discussion 

The nonperturbative part of the photon structure function has been of interest for a long 
time. At low it complicates efforts to check QCD experimentally and to extract the 
scale factor A. These problems can be dealt with in part by studying the photon structure 
function for higher moments where the prediction of QCD is stable against the higher order 



11 



corrections and where the nonperturbative part is small. (Of course this depends at least 
in part on the common assumption that the NP part vanishes quickly as the scaling variable 
X ^ 1.) The agreement of PQCD with experiment and a value A ~ 200 MeV is striking 
0. Nevertheless, it is still of interest to know to what extent the NP part of the structure 
function is really under control at low . 

We have argued here that one should explore the dependence of the virtual photon 
structure function on target . Then one can determine to an adequate extent just how 
important the NP part is. We think that this is clearly possible, given our results for the 
several NP models which we have studied. 

This has a number of consequences. First, by learning the actual size of -F/vp, we can 
improve checks on the PQCD prediction and extract a more reliable value for A. Secondly, 
the NP part and its relationship to the conventional PQCD prediction is of interest by 
itself. We know that checks of PQCD for high moments are satisfactory. However, we also 
know that we have no prediction for the n = 1 moment at = 0. It is a surprising and 
interesting fact that the virtual photon structure function has a calculable n = 2 moment. 
If we want to understand the relationship of perturbative and nonperturbative QCD, this 
problem needs to be explored experimentally. Equivalently, we need to understand the 
behavior of the structure function as x — > for both real and virtual photon targets. The 
existence of a finite or logarithmically rising cross section for the center of mass energy 
ly — > oo at fixed P^ ^ is a presumably nonperturbative phenomenon. This is related 
in some yet unknown way to the relation of perturbative and nonperturbative contributions 
to the structure function for small x. We think that this is an important issue. It is as yet 
unexplored by experiment. 

We can illustrate one unresolved issue by returning to our bland assumption that the 
NP part of the structure function can be obtained from the vector meson dominance model. 
While we have used this assumption to show that the NP part can be studied, we do question 
it. The P^ dependence can be studied by means of a dispersion relation in the target 



mass = fi , with threshold at /i = 4m^ 0, |Ty]. The NP part of the structure function 
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does have a contribution from small values of the mass parameter in this dispersion relation. 
In particular, there is a double pole at /i^ = corresponding to the simple p dominance 
model which we have studied. However, the NP terms in the photon structure function 
which are inverse powers of <^ A are not in fact so simple. We expect that the dominant 



terms behave as 10 



^2,Npi^J 1 — _|_ P2 ^2 ~'~ _|_ PL\ ' ^"^^^ 



In writing this, we have chosen to turn an expansion in powers of 1/P^ into one which is well 
behaved as 0. This represents an assumption about how higher order terms behave 

which we think is reasonable (but not proven). Only one of these terms can be thought of 
as corresponding to the vector meson double pole term, if mi nip. The other might have 
physically an interference term between a vector meson pole (for 1712 — > m^) and another 
term with no pole in /z^. We can see no reason why the second term cannot be important. 
We do not know of theoretical arguments which could establish its physical significance or 
size. 

We believe that these considerations make it clear why the dependence of the photon 
structure function is worth careful experimental study. We hope that these experiments 
will be carried out. As yet, little has been done in this promising area at the interface of 
perturbative and nonperturbative quantum chromodynamics. 
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Figures 



1. This figure set shows moments of the NP part of the photon structure function at 
p2 ^ 0, with Q"^ ranges from Q\ ^ 7i 45 GeV^ The moments for n = 2, 3, 4 
correspond to plots 1(a), 1(b), 1(c). Each curve contains NP behavior of the photon 
structure function for the three fits A, B, C mentioned in the text. 

2. These plots show the total real photon structure functions for the three fits. 2(a) and 
2(b) correspond to n = 4, 6 moments respectively. Each plot shows the corresponding 
total moment for the PQCD and NP structure functions superimposed. Dashed curves 
are PQCD parts and solid curves are total moments for the examples A, B, C. 

3. The figures 3(a), 3(b), and 3(c) correspond to n = 2, n = 4, and n = 6 moments of 
the virtual photon structure function respectively. There are also dependence in 
the figures, n — 1 case is the most interesting one since we can not see this moment in 
the real photon structure function. Dashed curves correspond to PQCD parts of the 
structure function and solid curves to total moments for the examples. A, B, C. 
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